Berry Phases with Real Hamiltonians With and Without a Many-body System as a 

Background 
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We present both the gauge theoretic description and the numerical calculations of the Berry phases 
with the real eigenstates, involving one with a many-body system as a background and the other 
with no such background. We demonstrate that for the former the sign of the Berry phase factor 
for a spin i particle (hole) coupled to a slow subsystem (phonon) depends on both the strength of 
electron correlations and the characteristics of the closed paths, unlike the cases for the latter. 
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The real wave functions in association with real Hamil- 
tonians are frequently encountered in numerous physi- 
cal problems, not to speak of condensed matter physics. 
There exists still an unresolved problem in defining the 
Berry phase concerned with the real eigenfunctions in 
that the gauge potential cannot be defined. The Berry 
phase arises from a [/(I) gauge potential as a result 
of the adiabatic transition involving the single- valued, 
complex, non-degenerate eigenstates of a fast subsystem 
coupled to a slow subsystem. In the following we briefly 
draw an attention to the main issue of the present work. 

In a system described by the time-dependent Hamilto- 
nian i/(R(t)) through a slowly varying parameter R, 



i/(R(t))|n;R(t)) =£;„(R(i))|n;R(t)) , 



(1) 



the total phase change of the state round a closed 
loop for a time period of T is given by 

|V(r)) =exp[z7„(C)]exp|-J^ dtE^{^{t))^ |^(0)) , 

where the Berry phase 7n(C) in the first factor is given 

by 



ln{C) 



I ^(n;R|V„|n;R) ■ dR , 



(2) 



and \ip{0)) — \n; R(0)) at t ~ 0. To obtain a non-zero 
value of 7„, the above eigenstate |n;R(i)) needs to be 
single- valued, complex and non-degenerate. Thus one 
cannot define the gauge potential 



A(R) = i{n;R\Vj,\n;K) 



(3) 



from the use of the real eigenfunctions. 

In order to allow for the case of multi-valuedness of the 
complex eigenstates \n; R) Berry derived the following 
expression of the 'magnetic field' from (ph. 



B(R) 



-Im 



t;R|(V^g)|m;R) x (m; R|(V^g)|n; R) 
[E,n{-R) - E^{R)Y 



(4) 



To avoid a confusion found in the literature, once again 
we stress that |n; R) above is the complex (but not real) 
and multi-valued eigenstate. From this expression one 
finds the existence of 'magnetic monopole(s)' correspond- 
ing to the degenerate eigenstates |m;R*) at singular 
point (s), R = R* in the parameter space R. 

The objective of the present study is two-fold; one is 
to rigorously discuss, in terms of the gauge potential, 
the Berry phase with the multi-valued real eigenfunc- 
tions, and the other, to newly examine the variation of 
the Berry phase with the strength of electron correlations 
for a coupled system made of a hole, the fast subsystem, 
and a phonon, the slow subsystem. In addition, we dis- 
cuss the two distinctively different cases from which the 
Berry phase arises; one from a many-body background 
system and the other from no such background. 

As mentioned above, the magnetic field B is not de- 
finable with the use of the multi-valued real eigenstates 
|n; R(i)), since the gauge potential A(R) vanishes. How- 
ever, by a proper gauge transformation it is possible to 
define the Berry phase as will be discussed below. The 
gauge invariance of the equation of motion (|l|) necessi- 
tates the local gauge transformation of the form 

|n;R)'=e'^"W|n;R) (5) 
The Berry phase 7n(C) is then 

In — i 'j>\n; RIVhI^; R)'- dR 

.R(T) 

= - / dA„(R) = -[A„(R(T)) - A„(R(0))] . (6) 



We now compute the magnetic field B(R), 

B(R) = V„ X A = X i\n\ RlV^jn; R)' 

= -Im ^XVfin-;R|TO;R)'xXm;R|VH|n;R)' (7) 



with m ^ n since "(Vj^ri; R n; R)' x '(n; R|Vj, |n; R)' is a 
real vector. We find from (|) 

\m-R\Va\7i-Ry= ^ " ^' , m^n. (8) 
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The insertion of (|) into (@) and a refinement of its result 
leads to 



B(R) = 



'Im 



;R|(V„i7)|m;R) x (m; R|(V«i7)|7i; R) 
[S„(R)-i;„(R)]2 



(9) 



Here both |m; R) and |n; R) are now the real eigenstates 
unlike the ones in Expression (^. It is important to 
realize from (|) and (|) that A(R) 7^ but B(R) = at 
other than the degenerate (singular) point(s) R = R*. 
Thus the Berry phase for the real eigenstates is realized 
as a type of the Aharonov-Bohm phase. 

Let us first examine a one-body problem involving no 
background. For a particle with spin in the constant mag- 
nitude of a magnetic field with slowly varying orientation 
constrained in the 2-D plane, the Hamiltonian becomes 

i?3 i?i 

matrix. Choosing A„(0) = nkO with k being a non-zero 
integer we write the gauge transformation, 



real; = cr ■ R 



with (T, the Pauli spin 



real eigenstates of |ri; 9) — 



In; 61)'= exp(m/c6i)|7i;6') . (10) 

The insertion of ([lo| ) into (|^) leads to the Berry phase 
of 7„ = —2'Knk. The eigenstate of an electron cou- 
pled to the magnetic field B confined in the 2-D x-z 
plane is real and double- valued. For the double- valued 

cos 0/2 \ . 1 
• /I /r, With n = as 
smt'/2 J ^ 

a spin component along the magnetic field, the Berry 
phase factor round the closed path which encircles k 
(many) 'solenoids' or a single solenoid of strength k is 
then 7n(C) = — /ctt. The Berry phase is thus equivalent 
to the flux in the 'magnetic solenoid' with 'strength' —n 
which vertically pierces through the 2-D plane if there ex- 
ists only one singular point in the plane. In the present 
case of a particle with no background, the Berry phase 
is independent of the choice of a path as long as the 
closed path encircles the solenoids. For the single unit 
of the solenoid, i.e., fc = 1, the Berry phase factor is 
7„(C) = — TT as expected. This is simply an Aharonov- 
Bohm phase type through the closed path which encloses 
the single solenoid fc = 1. It is now realized that the 
geometric phase (Berry phase) with the real eigenstates 
is composite in nature as it is determined by both the 
intrinsic property, i.e., the spin of the particle and the 
magnetic flux of the solenoid. Later we will examine 
the intrinsic property (spin) of a particle (hole) coupled 
to a phonon in the many-body background of antifer- 
romagnetic spin correlations. It is now clear that for 
the real eigenstates |n; R) the local gauge transformation 
n;R)'= e*'''"'-^^ |n; R) is essential for realizing the pres- 
ence of the gauge potential and the 'magnetic solenoid' 
as in the case of the Aharonov-Bohm phase. 

Let S* be a simply connected surface bounded by a 
closed loop in a parameter space R. If the phase change 



of the eigenstates \n\ R) occurs when transported adia- 
batically round a closed loop on S there must be at least 
one singular point where |n;R} is discontinuous, due to 
the intersection of potential energy surfaces . In the 
following we show an interesting case of the Aharonov- 
Bohm type Berry phase with the choice of a many-body 
system as a background, in which there can be any num- 
ber k of singular points through which the 'magnetic 
solenoids' pass vertically through the 2-D plane. 

By considering the Holstein-type electron (hole)- 
phonon coupling iQ-^ we use the effective single-band 
Holstein-Hubbard model Hamiltonian for a two- 
dimensional square lattice. 



nil 



<ij>,<T 



Y,n,R, + ^Y.i(^Tf + {^m (11) 



with Tlia- 



Here the real Hamiltonian i?(R) de- 



pends on the time-dependent external parameter R cor- 
responding to local lattice distortions. (cia) is the 
creation (annihilation) operator of an electron with spin 
a at site i. t represents the electron hopping strength; J7, 
the electron-correlation strength; 5, the electron (hole)- 
phonon coupling constant and K, the spring constant. 
The local Holstein distortion Ri for an in-plane breathing 
mode at the lattice site i — (ixjiy) is defined by 



A? 



-(1.0) 



-(04) 



(12) 



where Af (A^) is the displacement of oxygen from equi- 
librium along the positive x (y) direction at the lat- 
tice site i in the Cu02 plane and the other two, along 
the negative x (y) direction as shown in Fig. 1. Unlike 
the Holstein-i J model Hamiltonian [Tl|-[r^] the Holstein- 
Hubbard model Hamiltonian jl^ above has an advantage 
of investigating the geometric phase (Berry phase) vary- 
ing with the strength of correlation in the background. 




FIG. 1. Cu04 unit in the CUO2 planar square lattice to 
show a local lattice distortion (breathing mode). The solid 
circle denotes copper and the open circle, oxygen. 
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To obtain the Berry phase for the real ground state 
wave function, we treat a 4 x 4 periodic lattice with a 
hole surrounded by the background of nearly half-filled 
antiferromagnetic spins. To describe the propagation of 
the local lattice distortions round a closed path, we write 



(13) 



with r being the diniensionless time lapse, < r < 1 be- 
tween the lattice sites; the lattice distortion (breathing) 
occurs at site A at t = and at site i? at r = 1. In this 
manner, a closed path in the parameter space R can be 
defined. 

The mean field (Hartree-Fock) approximation ||l^ of 
(11) is 

H = -gY.4i^^R,{T)-tY^4i>, 

i <ij> 
i 

+[/^[(5,(z))2 + {S+{^)){SA^)) -{n{i)f] , (14) 



where the two component spinor is 



*T ) and the 



spin operators are defined by Sz{i) — 5(c||C^| ~ ^li'^iO^ 
= S-{i) = c\^c^^, and n{i) = c\^c^^ ^ 4^c^. 

By minimizing the total energy with respect to the 
average spin densities {Sz{i)), {S+{i)) and {S-{i)), the 
linearized Hamiltonian (|lj) was self-consistently treated 
for the lattice of the periodic 4x4 unit cell. For various 
values of U and g between Ot and lOt we evaluated the 
Berry phase factors. Schiittler et al. employed the 
Holstein-i J model by choosing sufficiently strong hole- 
phonon coupling. Our computed results of the Berry 
phase factor in the region of strong correlation and large 
hole-phonon coupling agreed with their results for the 
case of one-hole tunneling with various paths shown in 
Fig. 2. 

For numerical evaluations of the Berry phase factor 
we used the path integral approach |^,|l3|; the adiabatic 
quantum transition involving the eigenstate |n; R) round 
a closed path is given by 



Tnn{C) = exp 



Er,{B.{t))dt 



(n;R(T)|n;R(0)) , 
(15) 



where the first factor is the dynamic phase factor and the 
last one, the geometric phase (Berry phase) factor given 
by 



(n;R(T)|n;R(0)) = 



N 



lim V'(n;R(tfc)|n;R(tfe_i)) = e 

N^oo ^ — ^ 



(16) 



Only for the single- valued complex eigenstates |n;R) 
the Berry phase 7„ can be obtained directly from the 
relations (||) and (^. For the real eigenstate with the 
Hamiltonian (|lj) we evaluated the adiabatic Berry phase 
factor ( |l^ ) by computing the ground state electronic wave 
functions at each time step t^. Lately we [Tsl ] found that 
in scattering processes the role of non-adiabatic interme- 
diate transitions may be important. In such cases, the 
expression ( p^ ) is invalid. 



(a) C 





(c) 




(-) 




(+) 



FIG. 2. Hole hopping paths; the solid circle represents the 
initial location of the hole. Each sign denotes the Berry phase 
factor, +1 or —1. 



fc=l 



In Fig. 1 a single CUO4 unit |ll| , |12| that appears in the 
Cu02 planar square lattice of high copper-oxide super- 
conductors is drawn to exhibit the local lattice distortions 
A^'^ involving a copper atom at the center and its four 
surrounding oxygen atoms. By choosing the strong cor- 
relation oi U = Qt and the electron-phonon coupling of 
g — 6t for various closed paths, the sign of the computed 
Berry phase factors are shown in Fig. 2. In Fig. 3 we 
computed each single particle state involving 7 electrons 
with spin-up and 8 electrons with spin-down respectively 
for the doping case of one hole with spin-up for the path 
shown in Fig. 2(b). Due to the limited space we show only 
the lowest unoccupied energy level denoted as a dotted 
line for the spin- up electrons. A hole coupled to the lo- 
cal lattice distortion (breathing mode) may move around 
with some possibility of making a closed loop in the two- 
dimensional square lattice. A midway point (denoted as 
I in Fig. 3) is where a symmetric lattice distortion oc- 
curs between two neighboring sites, say, A and B. At 
such midway points the energy gap between the HOMO 
(the highest occupied molecular orbital) and the LUMO 
(the lowest unoccupied molecular orbital) is found to be a 
minimum, indicating that there exist degenerate points in 
the lattice distortion parameter space {Ri}. The HOMO- 
LUMO gap becomes a maximum at sites A, B and C 
which correspond to the symmetry breaking points. 

In the case of relatively strong electron correlations 
(large U), the computed Berry phase factor was —1 for 
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all the triangular closed paths and +1 for the closed path 
of the square shown in Fig. 2. We found that there exists 
an odd number of degenerate points (singular points), 
namely 3 of them for the former and an even number of 
degenerate points, i.e., 4 for the latter. For the case of 
sufficiently weak correlations, say, U < it is quite in- 
teresting to note that the computed Berry phase factor 
is +1 even for the triangular closed paths, contrary to 
the case of strong correlations. Thus our findings are as 
follows; the hole behaves as a particle of spin for the 
background of weak correlations and as a spin ^ fermion 
for the background of relatively strong correlations. We 
further found that the Berry phase factor for the hole 
which behaves as a spin i fermion in the background of 
strong correlations is subject to change its sign depend- 
ing on the nature of closed paths, i.e., the number of 
degenerate points that can exist round the closed paths 
in the parameter space. 



> 1 
c 






FIG. 3. Energy levels corresponding to lattice distortion 
along the closed path ABC A in Fig. 2(b). The dotted line in- 
dicates the lowest unoccupied energy level (LUMO) and the 
thick solid line, the highest occupied energy level (HOMO). 
Midway points (denoted as t) represent symmetric lattice dis- 
tortions between neighboring sites i and j = A, B, C). 

The direct use of the multi- valued real eigenstates for 
defining the 'magnetic flux' or the Berry phase is not 
justifiable. Wc argued, on the two grounds, the validity 
of the local gauge transformation |n;R)'= e'"^"^^) |n; R) 
with |n; R), the real eigenstate; namely the gauge invari- 
ance of the equation of motion (1) and the realization 
of A„ as a 'Aharonov-Bohm' phase type. In short, using 
the local gauge transformation we were able to define the 
'magnetic flux' of a solenoid and thus the Berry phase 
of the Aharonov-Bohm phase type. We then explored 
the Berry phase for the adiabatic transitions involving a 
hole coupled to the lattice distortion (breathing mode) 
in the antiferromagnetic background of both weakly and 



strongly correlated electrons. 

The following findings are in order. The strength of 
electron correlations affects the sign of the Berry phase 
factor; round the simplest triangular closed paths the 
Berry phase factor of the hole is -1-1 for the case of weak 
correlation, indicating the hole state of spin 0, and —1 for 
the case of strong correlation, indicating the hole state 
of spin i . However the Berry phase factor for the spin ^ 
fermionic hole depends on the nature of the closed paths, 
e.g., —1 for the triangular closed paths and for the 
square closed path. This is because the sign of the Berry 
phase factor is determined by the number k of degen- 
erate points (singular points) at which the fluxes pierce 
through the 2-D plane, depending on the characteristics 
of closed paths in the parameter space. Thus this is quite 
different from the system with no background in which 
the sign of the Berry phase does not depend on the choice 
of closed paths, as was discussed earlier. 
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